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Natural Capacity of a System of Two Two-Level Atoms 
as a Quantum Information Channel 
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A system of two closely spaced atoms interacting through a vacuum electromagnetic field is con- 
sidered. It is demonstrated that radiative decay in such a system resulting from photon exchange 
gives rise to a definite amount of information related to interatomic communication. Joint distribu- 
tions of detection probabilities of atomic quanta and the corresponding amount of communication 
information are calculated. 

PACS numbers: 03.67.-a, 03.65. B, 32.80.-t, 42.50.-p 



I. INTRODUCTION 



^ ' Analysis of physical systems as potential sources of quantum information is becoming an urgent issue in the context 

of extensive research into the physical implementation of quantum computation techniques [1-3]. Such an analysis is 
QQ • inevitably associated with mathematical problems that differ from those encountered in the analysis of such systems 
as objects of conventional methods of physical experiments. One of the natural goals of the approach specified above 
is to reveal the possibilities of using specific physical mechanisms that would ensure information exchange between 
microscopic quantum systems employed as components of quantum data converters. In this respect, two-level atoms 
^ (TLAs) interacting through a vacuum electro-magnetic field can be considered as one of the fundamental systems of 

this type. Obviously, two two-level atoms (TTLAs) form an elementary system. Certain efforts have been already 
concentrated on the investigation of such a system [4-7]. The main problem encountered in the exact calculation of 
^^ ■ TTLA dynamics stems from the necessity to rigorously take into account relaxation processes simultaneously with 
f^ , reversible interactions. 

^^ ■ In this paper, we investigate the process of radiative decay in a TTLA system in its pure form in the absence 

0^ ' of an external field. We consider TTLA dynamics on a time scale t ^ Tr, i.e., for time intervals greater than the 
p^ . radiative decay time r^ for a single atom. From the viewpoint of information transmission in a TTLA system, it is 
^ i' of interest to understand how much information is produced and stored in a system upon the completion of radiative 
', , decay processes in TLAs constituting the system under study and which type of information we deal with in this case. 
^ ' Evidently, such a formulation of the problem has no meaning for atoms separated by a distance on the order of or 
Ci ' greater than the wavelength, i.e., for atoms in traps, which can be considered as one of the prototypes of a physical 
^ ' quantum processor [3]. In this case, all the dynamic variables of a TTLA system decay on the same time scale r^.. 
However, for atoms of the relevant variables continue to relax within much greater time intervals t ^ Tr- Then, the 
information that relates the initial state of a TTLA system to its final state is stored. Such a geometry of a TTLA 
system may be implemented not only for atoms in dense media, but is also typical of impurity atoms adsorbed on a 
substrate. According to the experimental data of [8] , impurity atoms under these conditions may preserve the discrete 
C^ structure of atomic states and may be considered as microscopic quantum systems for quantum data con- version in a 
quantum computer. The mechanism behind the information exchange described above is associated with a relaxation 
photon exchange between closely spaced TLAs initially prepared in independent states. If the exchange rate 7cx is 
close to the radiative decay rate 7 of TLAs, then a single-quantum diatomic state ipa = |1) ® |0) — |0) (g) |1), which 
decays with a rate equal to 7 — 7ox, exists within a time interval t ^ r^ = 7^^. It is this singlet (in terms of a system 
of two spins) state that plays an important role in systems implementing the methods of quantum cryptography [9, 
10]. The mechanism under consideration naturally governs quantum data exchange between TLAs. 

In this paper, we restrict our consideration to a system of two closely spaced TLAs with a geometry of dipole 
moments shown in Fig. 1. Obviously, as long as the application of such a system is not specified, its information 
efficiency remains uncertain. Nevertheless, it seems appropriate to define the information efficiency in terms of the 
natural definition of the amount of data, which will be specified in this paper with the use of a given procedure 
of quantum measurement, without discussing the generalization of the considered approach. We assume that the 
information criterion employed in our study is sufficient for comparing different mechanisms of information exchange 
in terms of the information efficiency of feasible schemes for quantum computations using atomic transitions. In this 
paper, the definition of the amount of information will employ the Shannon definition of communication information 
corresponding to the probability distribution of measured energy quanta for each atom, which is described by a 
standard expression presented in [11]. 
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FIG. 1. Geometry of a TTLA system and relaxation processes of natural radiative decay. 



II. RADIATIVE DECAY IN A SYSTEM OF TWO CLOSELY SPACED ATOMS 



The system under consideration (Fig. 1) consists of two TLAs undergoing spontaneous radiative decay of excited 
states through the interaction with an electro-magnetic vacuum, which is described in a standard manner as a thermal 
reservoir. The relaxation superoperator £r of radiative decay can be calculated for a system of closely spaced atoms in 
the same way as in the case of a single atom [12] , i.e. , with a standard formula of second-order perturbation theory in the 
Hamiltonian of the interaction of an atom with positive- and negative- frequency components S,^^ = — cIJjEq (Rfe,T) 
of (Rfc are the coordinates of atoms) produced by emission and absorption of photons of the vacuum electromagnetic 
field Eg. Representing superoperators with the use of a substitution symbol of an operator being transformed, 
we arrive at the following expression for the superoperator C^ for several atoms in the Heisenberg representation (by 
analogy with the case of a single atom [13] )Q 
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Here, C^k is the relaxation eigensuperoperator for an isolated fcth atom due to spontaneous decay and >Crfcm describe 
relaxation photon exchange between the fcth and mth atoms, 
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Here ct^ = P12, (Tf, = P21, and jki is the relevant transition rate, 
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where n is the vector along the direction of the emitted photon, d^ is the transverse component of the dipole moment, 
Wa is the frequency of atomic transitions, and Rfem are the vectors of interatomic distances. In the case of two TLAs, 
along with the natural rate of radiative decay 7 for each atom, for the angle 9 between parallel dipole moments d]^2 
and R12 equal to 7r/2, we obtain a new relaxation constant of photon exchange 712 — ^cx' 
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where the modulus of the parameter g is less than unity, ip = uJs_Ri2/ c is the phase delay of a signal between the 
atoms. For 6 = 7r/2 and i?i2 -^ 0, we find that g k, 1 — 'pP' jh -^ 1, whereas in the case of anti-parallel dipole moments 
d-^2 J "^6 have g —* —\. Obviously, these two cases are equivalent to each other from the information point of view. 



^In the Schrodinger representation with a properly introduced Hilbert space of atomic operators, the corresponding operator 
is described by the superoperator £"*" Hermitian-conjugate of £. 



With the above-specified assumptions, the noise LiouviUian C^ is invariant with respect to the permutation of 
atoms and, correspondingly, has no matrix elements that would couple symmetric and antisymmetric operators of 
an atomic system. Therefore, the representation basis is chosen in such a manner that e-k = {s\ ® s|) /v^, k — 

1, ...,4; efe+4 = (si(^) ® s2^^,) + s\^^^^ ® s^^^)) /A fc = 1, ...,6, and eu+w = {s]^u) ® ^'(fc) " ^!(fe) ® ^'(fc)) /^' ^ = 
1, ..., 6, where a pair of indices i{k) and j{k) describes all possible combinations of pairs of single-atom basis vectors 
si and s|: {i(/c)j(fc)} = (12,13,14,23,24,34). The vectors s\ and s| are chosen as a basis of single-atom operators 

//•\/2, (P22 — Ai)/V2, (A2 + -P2i)/V2, and i(Pi2 — Ai)/'/2, orthonormalized with respect to a scalar product 
{A,B) = TrA+_B. Thus, the first ten basis elements describe symmetric atomic variables, whereas the last six basis 
elements describe antisymmetric atomic variables. Since, in the considered geometry, a laser field generates only 
symmetric excitations, the calculation of laser excitation is reduced to a ten-dimensional problem.g However, this 
paper does not consider a procedure of the preparation of the initial state. We are interested only in the transformation 



of the initial state of the form pi = p 
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in the process of relaxation. 



Applying superoperator (nj) to the basis elements e^, we obtain the following matrix representation for the relaxation 
operator: 
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where stands for a matrix with zero elements that complements the remaining submatriccs up to a 16x16 matrix, 
and 
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describe relaxation matrices for symmetric (S) and antisymmetric (A) variables. Expression (5) shows that these 
variables decay independently of each other, and small eigenvalues A^ ^ 7 may exist only for the Lgs matrix. 



^For g = 1, the dimensionality of the problem can be reduced to 9, i.e., to the dimensionality of a three-level quantum system 
that can be obtained from the starting four-level system by the exclusion of antisymmetric states. Generally, excitation with 
a vacuum leads to a random violation of the symmetry of wave functions but does not change this symmetry at the level of 
averaged fluctuations expressed in terms of the density matrix. 



III. DETERMINATION OF THE AMOUNT OF QUANTUM INFORMATION 

The final state at the moment of time t is described by a density matrix 

pj^S+if^^^p'^), (6) 

where S^ is the superoperator Hermitian-conjugate of the superoperator S — exp (£r^), which describes the transfor- 
mation of operators of physical variables in the Heisenberg representation. Here, we take into account two pairs of 
groups of dynamic variables — variables Ai and A2 related to atoms in the initial state at the moment of time t — 
before the interaction and the same variables at a certain final moment of time t — t after the interaction, which 
determines transform S. From the viewpoint of data exchange, it is of special interest to consider correlations between 
atomic variables at different moments of time, and r, rather than correlations between variables Ai{t) and ^2(t) 
at the same moment of time, which are described by expression (|6|) . We can calculate the information both between 
atomic TLA variables Ai(0) and A2(t) and between TLA pairs Ai(0) , ^2(0) and Ai{t) ,A2{t). 

To define the corresponding amount of information, we introduce the most typical class of measurement procedures 
that require the calculation of the relevant probability distributions. Consider a standard procedure of measurements 
with a coincidence scheme, which is described by a superoperator of the form I?+ = tt^ • • • ttJ" 7rJ~ • • • n~ , where tt^. are 
some Hermitian-field conjugate pairs of operators, in particular, operators of creation and annihilation. Suppose that a 
set of operators 7r,^(fc), fc = 1, . . . , n is defined for each m such that the condition of completeness ^j, 7r~ (fc)7r+ (fc) = / 
is satisfied. In fact, it is of interest to consider sets of eigenprojectors corresponding to transitions between n eigenstates 
?/^fc, Ami^k = ^ki^k, in accordance with relations 7r+(fc) = \ipk+i) (V'fcL where n + 1 — > 1 is taken for the brevity of 
notation. Here, we choose Am = fim for a single atom, where hm describes an operator of the number of an energy 
state at different moments of time t„i. For TLAs, we have irf — |2) (1|, tt^ — |1) (2|, tt^ = tt^, and ^2 = tt^. Now, 
we can introduce the relevant joint probability distribution, 

P(fci,...,fc„) = (P+(fci,...,fc„)p), (7) 

where the superoperator I?+(fci, ..., A:„) — 7r,^(fc„) • • ■ tt^ (ki) Q tti (ki) ■ ■ •7r~(fc„) introduces the corresponding measure- 
ment procedure, and averaging is performed over both TTLA and reservoir variables with allowance for the temporal 
evolution of the system under study within time intervals Tm+i — r^ between the moments of time r^ corresponding 
to the observation of detectable described by the set tt^. A formula 

£{ki,...,k„) =I?(fci,...,fc„)i 

relates the superoperator I?+(fci, ..., fc„) to a positive probability operator measure (POM) ||l[, which is conventionally 
employed in the quantum theory of information, i.e., a nonorthogonal expansion of unity £{dX), J £{dX) = /, which 
was introduced earlier in the theory of optimal quantum solutions (measurements) for the description of the procedure 
of a quantum measurement [14]. 

Here, we will use a set tt^ that corresponds to the populations of atomic levels at the moment of time ri = and 
an analogous set ttj corresponding to the moment of time T2 such that the inequalities 

7r2 > 1 > (1 - 5)7^2 , (8) 

are satisfied. In other words, we assume that the distance between the atoms is sufficiently small, so that the condition 
1 — g <C 1 is met. Then, if inequalities (^ are satisfied for a transient superoperator S — exp(£i.T2) we can employ 
the following approximation: 

S ^ lim exp(/:r |g=i T2) . (9) 

Under these conditions, the relaxation process brings a TTLA system into a stationary state, where the decay rate 
is equal to the rate of photon exchange, 57 = 7. In accordance with (8), this approximation is applicable on a time 
scale not greater than 

r=[{l-9h] - — 1^— j 7 , 
where A is the radiation wavelength for transitions in TLAs. 



Thus, we derive the following averaging rule for an operator of relaxation dynamics S with n = 2 in terms of 
the Heisenberg representation in the case when quantum numbers fci of the first atom are measured at i = and 
analogous numbers /c2 of the second atom are measured at t — t: 



Pi2{ki,k2) = TrpiT^ {ki)S [n^ (fc2)7r^(fc2)] T^'ti^i) ■ 



(10) 



With n — 4, when the quantum numbers ki and ^2 of the first and second atoms are measured at the moment t = 
and the numbers k^ and k^ of the same atoms are measured at i = t, we have 



P(i2){34:)(ki,k2,k3,k4) = TrpTT^ (^1)^2 (^2)5 [tt^ {k'i)^^ {ki)%^ {ki)%^ {k^)] tt} {k2)Tr^ (ki) . 
The amounts of information corresponding to these distributions are written as 

in the case of separate atoms and 

^(12)(34) = -ff(-P(12)) + H{P(34)) - ^(-P(12)(34)) 
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for pairs of atoms. Here, we employed a functional of the entropy H(w) = — ^ w(fc) log2 w(fc) and single-moment 

k 

probability distributions 



Pl(fc)=^Pl2(fc,fc2), P2(fc)=^Fl2(fcl,A:), 
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IV. RESULTS OF CALCULATIONS 

In the absence of relaxation, i.e., with 5 = 1, the final state m) coincides with the initial state, 7r2 = ttj^ , and joint 
distributions ( pl]| ) and ( fiA\ ) can be represented as products of two independent distributions of the form Pi{ki)P2{k2) 
and Pi{ki)P2{k3)S{ki — k2)S{k3 — k^j, respectively, where Pi{k) = Tr7r^(fc)/5^ )7rj~(fc). Obviously, the corresponding 
amounts of information described by ( |l2|) and (|l^) are equal to zero in this case, i.e., the establishment of information 
occurs only in the process of relaxation. 

Since we consider information in a quasi-stationary state, we have to calculate the transient superoperator (|9|) 
neglecting the decay of a single-quantum pure state ipa,, which occurs with a rate equal to (1 — g)j. To perform such 
a calculation, we use the following representation: 

S= \im\{£,-\iy\ 

where / describes the identity superoperator represented by an identity matrix of the corresponding dimensionality, 
which is equal to 16 in the case of TTLAs. The resulting matrix representation S with g = 1 is written as 
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For (7 < 1, we derive an expression 
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which takes into account the decay of all the atomic states except for the ground one. The density matrix of the 
stationary state in this case is written as 



p(i+2)(i) = 5+p(i+2) 



1 
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Such a density matrix corresponds to a quantum-free state '0o = |0) (E) |0), i.e., the ground state of each atom, and, 
consequently, contains no information concerning the initial state of atoms. By contrast, the stationary density matrix 
corresponding to expression (14) is written as 



/5(i+2)(f)=5+p(i+2' 
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where 
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Such a density matrix corresponds to a mixed state of the form 
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In this case, the quantity A represents the probability that a TTLA system resides in the singlet state i/^a upon the 
completion of fast relaxation processes. In fact, off-diagonal elements of the density matrix with indices 12 and 21 are 
not involved in the matrix corresponding to the joint probability distribution of the quanta of the first and second 
atoms: 
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FIG. 2. The amounts of information (a) between the initial state of the first atom and the final state of the second atom 
and (b) between the initial and final states of both atoms as functions of the mean populations of excited levels of the initial 
atomic states. 



For g = —1, we obtain exactly the same distribution. The corresponding amount of information (y_2|) is presented in 
Fig. 2a. The maximum of the information obtained between the final state of the second atom and the initial state 
of the first atom is achieved with a population n2 = or n2 = 1 for the second atom and ?t.2 = — 1 + V37/4 = 0.521 
or, respectively, ^2 = 2 — v37/4 = 0.479 for the first atom. The maximum amount of information is given by 
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and is approximately equal to 0.14 bit. 



Joint two-atom two-moment distribution ( ]11| ) is represented by a matrix 
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The corresponding amount of information ( |l3|) in both atoms at the moment t = r as compared with the initial state 
is presented in Fig. ||b. The maximum of this quantity is ^(™2V34') ~ 0.322 bit. 

V. CONCLUSIONS 

The results of calculations presented above demonstrate that there exists a natural mechanism of information 
exchange between two closely spaced two-level atoms. If these atoms are initially prepared in statistically independent 
states, then, after a lapse of time r such that (1 — g)^T <C 1 ^ 7r, where 7 is the decay rate of single-atom excitations 
and g « 1 is the efficiency of interatomic photon exchange with respect to single-atom decay, the procedure of counting 
the quanta of atomic excitations reveals the presence of nonzero information between the population of the second 
atom at the moment of time t and the initial population of the first atom. The maximum amount of information 
-^max — 0.14 bit corresponds to the upper or lower initial state of the second atom and the mean population of the first 
atom close to 0.5. An analogous amount of information in both atoms is equal to /max = 0.322 bit. Thus, a certain 
amount of information is produced without special efforts due to the specific features of spontaneous radiative decay 
in a system of closely spaced atoms. 

Taking into account not only the relaxation but also resonant electrostatic dipole-dipole interaction, which is 
described by a Hamiltonian of the form Tic — Qcipt ^^. + ^\^2^^ 0116 can easily verify by direct calculations 
that the superoperator of the transient distribution dij) remains unchanged. The same conclusion follows from the 
consideration of the symmetry of the system with allowance for the invariance of 1-ic with respect to the permutation 
of atoms. Thus, electrostatic dipole-dipole interaction has no influence on relaxation data exchange as long as this 
interaction remains too weak to change the relaxation operator. 

The results obtained above cannot be directly applied to the analysis of atoms on a surface because of large 
nonradiative dephasing rates, which gives rise to the main difficulty encountered when attempts are made to employ 
adsorbed atoms for the implementation of quantum computations [3]. However, we should note that, for closely spaced 
atoms, analysis of dephasing effects should take into account the specific features of relaxation processes associated 
with the interaction of atoms through dephasing excitations. The mechanism behind such an interaction should be, 
in a certain respect, similar to the mechanism considered in this paper. 
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